Spatially inhomogeneous phase evolution of a two-component Bose-Einstein 

condensate 
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We investigate the spatially dependent relative phase evolution of an elongated two-component 
Bose-Einstein condensate. The pseudospin-i system is comprised of the \F = l,mp = —1} and 
\F — 2,m,F = +1) hyperfine ground states of 87 Rb , which we magnetically trap and interrogate 
with radio-frequency and microwave fields. We probe the relative phase evolution with Ramsey 
interferometry and observe a temporal decay of the interferometric contrast well described by a 
mean-field formalism. Inhomogeneity of the collective relative phase dominates the loss of inter- 
ferometric contrast, rather than decoherence or phase diffusion. We demonstrate a technique to 
simultaneously image each state, yielding subpercent variations in the measured relative number 
while preserving the spatial mode of each component. In addition, we propose a spatially sensitive 
interferometric technique to image the relative phase. 

PACS numbers: 03.75.Kk, 03.75.Dg, 03.75.Mn, 34.50.Cx, 42.30.Rx 



I. INTRODUCTION 

Knowledge of the phase of matter waves is crucially im- 
portant in studies of interferometry, entanglement, and 
precision measurement. Pertinent to all of these areas is 
the pseudospin-i condensate, which can be realized with 
trapped neutral atoms in a superposition of two hyperfine 
ground states. Most commonly, the magnetically trap- 
pable \F = l,m,F = —1) and \F = 2,mp = 1) states of 
87 Rbhave been used in experiments [1, 2, 3, 4, 5]. Early 
studies demonstrated spatial separation of the compo- 
nents [1] and interferometric detection of relative phase 
for a region where both components remained overlapped 
following strongly damped center of mass motion [2]. 
Also, the effect of phase winding throughout a two- 
component condensate was studied for continuous elec- 
tromagnetic coupling [3]. Following this, images of spin 
excitations in an ultracold uncondensed gas [6] enabled 
the study of spin domain growth for mixtures of con- 
densed and uncondensed atoms [4]. The interference 
between two vortex lattices comprised of each compo- 
nent has also been examined [7]. More recently, long- 
lived ringlike excitations of the binary condensate sys- 
tem have been observed [5] . This kind of dynamical in- 
stability is accompanied by spatially dependent relative 
phase dynamics, which we investigate here. In particu- 
lar, we consider the temporal decay of the interference 
signal obtained with a Ramsey-like measurement of the 
pseudospin-i condensate. 

The mechanism of phase diffusion for a two-component 
quantum gas has been recently studied [8], whereby the 
evolution of a coherent spin state results in the decay of 
Ramsey visibility [9]. For the close inter- and intrastate 
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interaction strengths in our system, phase diffusion is 
negligible [9] . Rather, we consider mean-field driven spa- 
tial inhomogeneities of the relative phase, which also act 
to decrease the interferometric contrast even without sig- 
nificant spatial separation between the components. This 
is relevant in the context of proposals to squeeze the 
macroscopic pseudospin in two-component quantum de- 
generate gases [10, 11, 12, 13], which are based on using 
the mean- field interaction as a source of entanglement, 
and development of a trapped atomic clock using an atom 
chip [14, 15]. 

This paper is organized as follows. In Sec. II, our 
experimental procedure is described up to the point of 
initialization of a two-component Bose-Einstein conden- 
sate (BEC). In Sec. Ill, we report on the spatiotemporal 
relative phase dynamics of the two-component conden- 
sate using Ramsey interferometry. Simulations of cou- 
pled Gross Pitacvskii equations including atomic loss and 
electromagnetic driving terms yield striking agreement 
with observed matter wave interference. In Sec. IV, we 
present the demonstration of a simultaneous state selec- 
tive imaging technique. The method enables improved 
measurement of the longitudinal spin projection in a sin- 
gle experimental run, when classical fluctuations in rela- 
tive and total atom numbers exist between different re- 
alizations of an experiment. Finally, in Sec. V, we pro- 
pose a combination of the dual state imaging method and 
Ramsey interferometry to directly image the spatial in- 
homogeneity and quantum fluctuations of relative phase. 



II. REALIZATION OF TWO-COMPONENT BEC 

A detailed description of our apparatus, including ex- 
periments performed with a \F = 2,mp = 2) condensate 
on a perpendicularly magnetized film atom chip, has been 
described elsewhere [16, 17, 18]. In brief, we use an 
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FIG. 1: (a) Hyperfine ground states of 87 Rb . Two Zeeman 
levels are coupled by a two-photon microwave radio-frequency 
field (see text for description), (b) Two-photon Rabi oscilla- 
tions of a pseduo-spin- 1 BEC as measured by simultaneous 
detection of the populations of |1) and |2) (Sec. IV). 



atom chip with a machined Ag foil structure that allows 
currents to be passed in U- and Z-shape paths for sur- 
face magneto-optical trapping and Ioffc-Pritchard mag- 
netic trapping [19]. We begin by preparing a condensate 
in the \F — l,m,F = — 1) state, hereafter referred to as 
state |1). This is achieved in a manner similar to previ- 
ous work with two main differences. First, we optically 
pump the atoms into |1) prior to magnetic trapping us- 
ing a 1 ms duration pulse of <j~ polarized light, tuned 
to the D 2 (F = 2 — > F' = 2) transition, in the ab- 
sence of repumping (F = 1 — * F' = 2) radiation. To 
ensure the purity of |1) during the magnetic trapping 
stage, a 2 ms pulse of optical pumping light illuminates 
the trapped cloud, completely removing residual mag- 
netically trapped atoms in the F = 2 level. Second, the 
BEC is imaged after magnetic trapping via optical ab- 
sorption using a 100 /is pulse of a + light, tuned to the 
D 2 (F = 2 — > F 1 = 3) transition. This is immediately 
preceded by a 1 ms pulse of repumping light that trans- 
fers all the atoms into the F = 2 manifold for imaging, 
while the short repumping- imaging delay ensures image 
blurring (from the recoil of a single repumping photon 
by each atom) is minimized. A charge-coupled device 
(CCD) camera records the absorption image using an 
achromat doublet lens, with a resolution of 7.5 /im/pixel. 
Using this sequence, a pure |1) BEC of 2 x 10 5 atoms is 
routinely created in a cycle time of 40 s. 

Using a |1) condensate, the radial and axial trap 
frequencies of the magnetic trap were measured to be 
f p = 97.6(2) Hz and f z = 11.96(2) Hz. We observe a |1) 
condensate lifetime of 1.7 s, well below the observed mag- 
netic trap lifetime for a |1) thermal cloud (30 s) and the 
expected three-body collisional-loss- dominated lifetime 
(17s). Since the total atom number, BEC plus ther- 
mal cloud, is conserved, we attribute the observed life- 
time to be limited by heating due to technical noise. 
The rate of this heating was measured to be 12(1) nK/s 
which is small compared to the transition temperature 
of 120 nK, yet ultimately of consequence. We prepare 
a two-component condensate of 1.5 x 10 5 atoms using 



a two-photon microwave radio-frequency field to couple 
|1) to \F = 2,toj? = 1), hereafter referred to as state |2). 
This is schematically presented in Fig. 1(a) and was ini- 
tially demonstrated and described in Ref. 20. We use a 
magnetic trap with the potential minimum at 3.23 G (for 
which the relative Zccman shift is independent of mag- 
netic field to first order), as determined by microwave 
spectroscopy of a trapped condensate. The ~ 6.8 GHz 
microwave radiation is derived from a signal generator 
(Agilent E8257D), pulsed using a fast ns switch (Agilent), 
then amplified (MA. Ltd 10 W) and transmitted to the 
BEC using a helical antenna (gain 19 dB) which resides in 
air 10 cm from the condensate. We simultaneously drive 
a radio-frequency magnetic field at ~ 2.01 MHz with the 
atom chip wires used for evaporative cooling, albeit with 
an independent signal generator (SRS DS345) and an 
amplifier (MA. Ltd 2W). These fields couple |1) and |2) 
via the intermediate state \F = 2, mp = 0) with a detun- 
ing of 250 kHz. This detuning ensures that no significant 
population is transferred to \F = 2,mj? = 0), a magnet- 
ically untrapped state, such that a two-level formalism 
adequately describes the multi-component dynamics of 
the system. We measure a two-photon Rabi frequency 
of 416 Hz [Fig. 1(b)] which is much faster than the char- 
acteristic time scale for the condensate to change shape. 
The high quality of these data is a result of the simulta- 
neous state detection method described in Sec. IV. 



III. RAMSEY INTERFEROMETRY USING 
TWO-COMPONENT BEC 



Here we use Ramsey interferometry to study the rela- 
tive phase evolution of a two-component condensate. The 
non-equilibrium dynamics observed in this system [5] are 
accompanied by spatially dependent dynamics of the rel- 
ative phase. We describe the evolution of the system us- 
ing a pscudospinor formalism, whereby the internal and 
external states of the condensate are represented by a 
two-component order parameter 
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where n\ and n 2 are the atomic densities of each state in 
the condensate. We begin with a |1) condensate in the 
ground state of the combined mean-field and external 
potentials, with density no(r). This corresponds to the 
pseudospinor representation 
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Application of the first 7r/2 pulse (of length t^/ 2 ) prepares 
the two-component superposition [21] 
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This is no longer the ground state of the two-component 
system, as the mean-field interaction between component 
|1) with component |2), and component |2) with itself is 
different to that of component |1) alone. This is due 
to a slight difference in the s-wave scattering lengths 
an = 100.40 ao, a.22 = 95.00 ao, and an = 97.66 ao, 
where ao is the Bohr radius [5] . Through the state inter- 
conversion, we modify the mean-field energy of the sys- 
tem by ~ 1% for our experimental parameters and a 
condensate with 1.5 x 10 5 atoms. This is enough to drive 
hundreds of milliseconds of weakly damped collective ex- 
citations and coherent relative phase evolution. The sys- 
tem is allowed to evolve for a time T, after which we 
observe the condensate with or without the application 
of a second 7r/2 pulse 
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where primes denote the densities and phases after the 
optional second pulse. We define a spatially dependent 
relative phase as 



(r) = 4>2{r) - c/)i(r) 



(5) 



In the frame rotating at the effective frequency of the 
two-photon coupling, this phase evolves to 



(r) = AT + mf (r) + ^ 



(6) 



after evolution time T, where A is the detuning of the 
two-photon field from the atomic resonance, </> m f is the 
spatially dependent phase whose evolution is driven by 
the mean field, and 5(f) l& anv phase shift intentionally ap- 
plied to the coupling field before t = t n / 2 + T. Varying A, 
T, or 8(j> and applying a second n/2 pulse yields a Ramsey 
interference signal a, which we define as the expectation 
value of the longitudinal psuedo-spin projection, scaled 
to vary as the normalized population difference 



(Nh)- 1 (S Z ) = (N^-N[)/N , 



(7) 



where N! = J n^dr, (i = 1,2) and N is the total atom 
number. 

A typical Ramsey signal obtained by varying the evo- 
lution time T is shown in Fig. 2. We observe nonexpo- 
nential decay and chirp to Ramsey fringes in the time 
domain. These observations are consistent with the spa- 
tial evolution of a coherent two-component order param- 
eter, as described below. Decoherence does not signif- 
icantly contribute to the decay of Ramsey contrast we 
observe [22]. It is the spatial dependence of the mean- 
field driven phase that leads to a decrease in contrast of 
Ramsey fringes over time. It should be noted that this 
dephasing is unlike the inhomogeneous dephasing that 
occurs in statistical ensembles of uncondensed thermal 
atomic samples. Rather, it is the variation of the relative 
phase across a many-body wave function. 
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FIG. 2: The Ramsey interference signal of a pseudospin 87 Rb 
BEC where the detuning of the two-photon field from the 
atomic resonance A/2tt = 1Hz. The points are experimental 
measurements and the solid line is calculated using the CGPE 
theory, with an additional exponential decay due to decoher- 
ence [22]. The ~ 16 Hz oscillation, decay, and frequency chirp 
of a are predominantly due to spatial evolution of the relative 
phase. 



To simulate the evolution of density and relative 
phase, the dynamics throughout electromagnetic cou- 
pling pulses, and the Ramsey interferometry signal, we 
solve the coupled Gross-Pitaevskii equations (CGPE) 
for both components with decay terms corresponding 
inter- and intrastate many-body loss processes. In their 
most basic form (without loss or electromagnetic cou- 
pling terms), these equations were introduced in [23, 24]. 
The inclusion of electromagnetic coupling terms has ap- 
peared in [25] (multicomponcnt condensates) and [26, 27, 
28, 29, 30, 31] (specifically pseudospin-i condensates). 
Addition of nonlinear terms to describe many-body loss 
processes was used in [32] and applied to a binary con- 
densate system of 87 Rb in [5] . Combining all of the above 
features, we arrive at 
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where m is the mass of 87 Rb , Vi is the magnetic 
trapping potential experienced by component i, and 
cjij = ATth 2 a.ij /m are the mean- field interaction param- 
eters. The dominant two- and three-body loss processes 
are described by the terms Ti = (7mn 2 + 7i 2 n 2 )/2 and 



T 2 = (7 22 n 2 +7i 2 ni)/2, with loss rates ji, 
and described in Rcf. 5. 
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FIG. 3: Simulated in-trap column densities J m dx and rela- 
tive phase 4>(p = 0, z) [see Eqs. (4) and (5)] for various evo- 
lution times of the initial state in Eq. (3). The initially uni- 
form phase exhibits inhomogeneous evolution, accompanied 
by density excitations in the axial direction (z). An offset of 
14 7rrads~ has been added to the relative phase for clarity. 



To solve Eqs. (8) numerically, we exploit the 
cylindrical symmetry of our experimental geometry, 
and assume that the wave function takes the form 
^(f) — ► e lmv \E' m (p, z), where m is an integer. Since the 
ground state of the system has m = 0, we exclusively 
solve for this case and assume the time evolution does not 
break this symmetry We use the discrete Hankcl-Fouricr 
transform to solve the CGPE and have adapted the pro- 
cedure used in [33] for multicomponcnt systems. Essen- 
tially, the technique facilitates the solution of second- 
order partial differential equations in cylindrical coordi- 
nates using spectral methods, where the Fourier trans- 
form fails due to divergence of the Laplacian for zero 
radial momentum. Implementation of the discrete Han- 
kel transform is greatly simplified by sampling the fields 
at zeros of the m th order Besscl function of the first kind, 
rather than using a Cartesian grid. 

In Fig. 3, we show the simulated column densities 
of each component for various evolution times following 
preparation of the initial superposition in Eq. (2). The 
relative phase is seen to vary across the axial dimension 
of the condensate. This relative phase modulation even- 
tually becomes manifest as relative density modulation, 
such that the components spatially separate along the 
axial direction, consistent with the sign of 011022 — o>f 2 
that governs the miscibility of two-component conden- 
sates [23]. This is evident in Fig. 3 at 90 ms, where 
|1) predominantly resides further from the origin than 
|2). In the absence of loss, the moments J \z\rii df would 
continue to oscillate out of phase, with |2) occupying a 
larger volume than |1) at ~ 200 ms. However, atomic loss 



dampens the collective oscillations on this time scale. 

We present here measurements of the density of state 
|2) for various evolution times T, followed by a second 
7r/2 pulse. We observe striking agreement between these 
results and those of CGPE simulations. In Fig. 4, the 
column density J 712 dx is shown, where x is the imaging 
axis, a component of the radial coordinate (of tight con- 
finement) p. The only free parameter used in the sim- 
ulation is the detuning A of the two-photon field from 
the atomic transition. We adjust this parameter to op- 
timize the agreement (in this case A/2ir = 1Hz), con- 
sistent with a measurement of the detuning using a cold 
thermal cloud. In particular, the three distinct density 
maxima at 40 and 115 ms, corresponding to antinodes of 
the state |2) wave function, are represented accurately by 
the simulation. This prominent density modulation along 
the z axis (of weak confinement) reveals the nonuniform 
relative phase variation along this direction. This is en- 
abled by the second 7r/2 pulse, which converts informa- 
tion about the relative phase into (relative) density. In 
the absence of the second pulse, while density modula- 
tion still exists, it is of lower contrast (Fig. 3, 90 ms is 
the most structure we observe from simulations). 

Furthermore, the Ramsey contrast decays significantly 
faster than the components spatially separate. We em- 
phasize that uncommon spatial modes of the components 
are neither necessary nor sufficient for loss of interfero- 
metric visibility. For example, consider the ground state 
of the two-component system, whose density consists of 
a central feature of state |2), enclosed by two lobes of 
state |1). Of course, no inhomogeneous relative phase 
evolution occurs for the ground state and the mean-field 
limited Ramsey visibility is 80%. This remains true even 
in the presence of loss, which we have confirmed from 
simulations that show the condensate wave functions adi- 
abatically follow the time-dependent ground state. 

We attribute residual differences in the precise spa- 
tial structure of the measured and simulated macroscopic 
wave functions to the finite resolution of our imaging sys- 
tem, photon recoil blurring during the imaging pulse, im- 
perfect cylindrical symmetry of the trap, and interaction 
of the condensate with an undetectable yet significant 
thermal component for longer times. The repeatability 
of the observed interference fringes (both run-to-run and 
day-to-day) is indicative of the well-defined initial-state 
preparation and phase evolution in this system. 



IV. DUAL STATE IMAGING 

We present a new imaging technique to simultaneously 
detect states |1) and |2) with spatial resolution. This 
represents the first unambiguous state selective measure- 
ment of a pseudospin-i condensate that preserves the 
spatial mode, pseudospin-^ condensates of 87 Rb suitable 
for intcrfcromctry are commonly comprised of states 
whose magnetic moments are similar, complicating spa- 
tial separation of the constituent states with a magnetic 
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FIG. 4: Column density of state |2) as a function of Ramsey interferometry time (inset for each image in milliseconds). 
Application of the second tt/2 pulse locally converts the relative phase into relative population. (Left) Experimental: single 
shot absorption images, taken after a 20 ms fall time, reveal the mean-field driven phase evolution. (Right) Theoretical: column 
density plots obtained by solving the coupled Gross-Pitaevskii equations show excellent agreement with experimental results. 
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FIG. 5: Transfer of |1) to the magnetic field insensitive state 
|2,0) is achieved by adiabatically sweeping the atomic reso- 
nance in the presence of a fixed frequency microwave pulse, 
(a) Spatial separation is performed by pulsing I z (see text for 
description) to induce a Stern-Gerlach force that accelerates 
1 2) away from the atom chip and |1). (b) Absorption imaging 
using the D2 (F = 2 — > F' — 3) transition yields a single 
shot image of a two-component condensate, |1) (above) and 
1 2) (below), where the |1) structure arises during a Ramsey 
interference experiment described in Sec. III. 



field gradient. Various methods exist to individually im- 
age either the F — 1 or F = 2 populations, using reso- 
nant light in combination with optical or microwave re- 
pumping pulses. An alternative method uses a radio- 
frequency sweep to distribute the |1) and |2) populations 
among all eight Zeeman sublevels which are then spa- 
tially resolved using Stern-Gerlach separation prior to 
imaging. The original |1) and |2) populations are in- 
ferred from estimates of how the adiabatic rapid passage 
distributes population among the sublevels [34]. Phase- 
contrast imaging can also been used to image the differ- 
ence in populations [3] and has enabled the observation 
of vortex dynamics [7, 35, 36]. 

Our dual state detection method consists of three 



stages. First, the magnetic trap is switched off by re- 
ducing the Z-wire current, I z , to zero in ~ 200 /is. The 
uniform bias field that completes the trapping potential 
is concurrently ramped off in 3 ms and a fixed frequency 
microwave field is applied. The Zeeman splitting of the 
magnetic sublevels decreases, facilitating rapid passage 
from |1) to \F = 2,mp = 0). This population transfer 
method is more robust than a resonant microwave 7r 
pulse, owing to the less stringent requirements of rapid 
passage on the stability of magnetic fields and microwave 
power. The same microwave field is used for the two- 
photon pulses, obviating the need for a second microwave 
frequency synthesizer and enabling detection immedi- 
ately after a two-photon pulse, as in Sec. III. For our ini- 
tial experiments demonstrating this technique, we trans- 
ferred 44(1)% of |1) population to \F = 2,m F = 0) (rapid 
bias field switching was desired at the expense of adia- 
baticity). At this point, both spin components of the 
condensate are spatially overlapped and fall freely under 
gravity. However those in the \F = 2, mj? = 0) state are 
largely insensitive to magnetic fields. 

The second step [Fig. 5(a)] produces spatial separation 
of the two components. This is achieved using a 30 A, 
1 ms I z pulse which induces a Stern-Gerlach force in the 
direction of gravity to displace |2) such that both compo- 
nents are separated by ~ 300 /im after 20 ms of free fall. 
Finally, an absorption image is taken [Fig. 5(b)] yielding 
the column density of |2) directly below |1). This image 
is from a Ramsey interferometry experiment described 
in Sec. Ill and highlights the ability of this dual state 
imaging method to preserve spatial information about 
the wave functions. This is due to the minuscule recoil 
velocity (5 x 10 -9 mm/s) associated with microwave pho- 
ton absorption thus eliminating the image blur associated 
with optical rcpumping techniques. 

While conservation of spatial information is crucial 
for dynamical studies of two-component BEC, the de- 
tection method should also determine the atom number 
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FIG. 6: (a) The detection efficiency 77 for dual state absorp- 
tion imaging as a function of imaging pulse duration, for three 
typical saturation parameters, s — 0.50 (solid), s = 0.10 
(dash-dotted) and s = 0.01 (dashed), (b) Simultaneous mea- 
surement of atoms in |1) and |2) after a two-photon coupling 
pulse showing high correlations for a range of detected BEC 
atom numbers. An imaging pulse of 100 ps duration with 
s = 0.1 was used. 



accurately. Wc investigate the effects of optical pump- 
ing during absorption imaging on r), the relative effi- 
ciency of detecting atoms originating in \F = 2, mj? = 0) 
to those starting in |2). The absorption of a <j + polar- 
ized probe beam was modeled using the optical Bloch 
equations [37], including the Zeeman splitting of the 
ground and excited states. We evaluated the absorp- 
tion by integrating the total excited-state population 
(F' = 1, 2, 3 ; m' F = —F',,.,F') during the imaging pulse 
for a given initial state. Results arc shown in Fig. 6(a) for 
different saturation parameters s = Jo/1.68mWcm -2 of 
the \F = 2, mp = 2) — > \F' = 3, m' F = 3) cycling transi- 
tion, with Iq the intensity of the probe beam. As the 
imaging pulse duration increases, optical pumping acts 
to populate the \F = 2, mp = 2) state and more time is 
spent driving the cycling transition. As such, the rel- 
ative detection efficiency asymptotes to unity and does 
so faster for higher imaging intensities. For our typical 
imaging conditions, we predict a relative detection effi- 
ciency of 77 = 0.98. 

The dual state imaging method minimizes fluctuations 
in measured relative number that arise when detecting 
individual populations in different experimental realiza- 
tions. This is due to common mode rejection of both 
probe-laser frequency noise and irreproducibility of total 
atom number between experimental runs. This is demon- 
strated in Fig. 6(b) which shows a strong correlation be- 
tween |1) and 1 2) populations for a wide range of detected 
total condensate atom number. The method allows the 
relative population to be measured with a subpercent 
standard deviation, as opposed to single-component mea- 
surements of the absolute population which can typically 
vary by ±10%. This method improves our measurement 
of phenomena such as Rabi oscillations [Fig. 1(b)] and 
Ramsey fringes (Fig. 2). Moreover, it will see impor- 
tant application in atom chip clock research, allowing 



normalization of total atom number fluctuations, a prac- 
tice used to enhance the accuracy of precision fountain 
atomic clocks [15, 38]. It also represents a significant step 
to achieving quantum-limited detection of the total spin 
projection in pseudospin-i condensates. 



V. PHASE RECONSTRUCTION USING 
RAMSEY INTERFEROMETRY 

In this system, imaging the relative phase amounts to 
resolving the transverse component of the local Bloch- 
vector corresponding to the pseudospinor in Eq. (1). 
For a sample of ultracold uncondensed 87 Rb atoms, this 
has been achieved in Ref. 6. There, the Ramsey inter- 
ference signal was binned in both space and time, re- 
covering a phase spectrogram that illustrated collective 
spin excitations driven by exchange scattering. This re- 
quires the repetition of many experimental cycles and 
gives an averaged relative phase with reduced tempo- 
ral resolution. We propose a technique to spatially re- 
solve the relative phase of the pseudospin-i condensate 
from two dual state absorption images and present the- 
oretical results that investigate the method's applicabil- 
ity and limitations. Performing an identical experiment 
with and without the presence of the second 7r/2 pulse 
permits four column densities h = J n dx to be imaged, 
where n € {ni, U2, n'i, n' 2 } and the densities are defined 
in Eqs. (4). The densities can be used to recover the 
relative phase <f>(r) using 



n' 2 (r) — n[(r) = 2y/ ni{r)n 2 (r) sin(</>(r)) 



(9) 



There is little variation of the relative phase along the 
radial direction, and the densities of both components 
have similar radial dependence. This is a consequence of 
the tight radial confinement, as the energy of the system 
is not sufficient for significant radial excitations to be 
induced. As such, the radially integrated densities, the 
line densities n = J n2irpdp = J ndy, are representative 
of the relative phase evolution, such that 

n' 2 (z) - n'^z) = 2^n 1 {z)n 2 (z)sm(4>(z)) (10) 

defines a phase 4>(z) that is a good estimate of the phase 
along the center of the condensate <f>{p = 0, z). We have 
confirmed this by evaluating Eq. (10) with line densities 
attained from simulations, as shown in Fig. 7. For 90 ms 
of evolution, where the range of relative phase is 1.9 tt 
across the condensate, the agreement is within 78 mrad. 
In more elongated geometries or where tunable miscibil- 
ity can be realized, radial excitations are further sup- 
pressed and the approximation improves. In these cases 
there is increased spin locking along the imaging axis and 
the integration of density inherent to absorption imaging 
corrupts the phase signal less. 

The requirements of this technique are more demand- 
ing than those of the experimental investigation pre- 
sented in Sec. III. For successful implementation of 
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FIG. 7: Simulation of interferometric phase retrieval method. 

(a) Two absorption images are captured for each evolution 
time, yielding the line densities immediately prior to, and 
after the second 7r/2 pulse: fix, fi2, and n[, n'2 respectively. 

(b) Point-by-point arithmetic is used [Eq. (10)] to extract the 
relative phase profile along the axial direction. 

Eq. 10 to point-by-point arithmetic on multiple absorp- 
tion images, we are currently pursuing improved detec- 
tion methods with higher signal-to-noise ratio, spatial 
resolution, and dynamic range. This technique may be 
useful for investigations of spin squeezing, in the context 
of recent proposals which suggest using the nonlinear in- 
teraction to reduce the spin projection noise below the 
standard quantum limit [10, 11, 12, 13]. It is interesting 
to consider the spatial dependence of the relative phase 
fluctuations of a macroscopic pseudospinor. This tech- 
nique may allow such quantum noise to be studied with 



spatial resolution. 



VI. CONCLUSION 

We have quantitatively studied the spatial evolution of 
the relative phase of a pseudospin-i 87 Rb BEC. Prepar- 
ing the system with a 7r/2 pulse yields a non-equilibrium 
state due to small differences in the inter- and intra-state 
scattering lengths an, 022, and a 12 . We demonstrate the 
resulting spatially inhomogeneous evolution of the col- 
lective relative phase. This acts to decrease the Ram- 
sey contrast before significant spatial separation of the 
components. Spatial dephasing of the many-body wave 
function dominates the effects of decoherence and phase 
diffusion on the interferometric contrast. A second 7r/2 
pulse maps the non-uniform phase distribution onto the 
population distribution of each spin state. Measuring 
the longitudinal pscudospin projection is enhanced by a 
new dual state imaging method that combines microwave 
adiabatic rapid passage and Stern-Gerlach spin filtering. 
This yields a single shot measurement of |1) and |2), en- 
abling the relative populations to be measured with sub- 
percent precision. Moreover, this novel method does not 
corrupt the spatial modes of either component. We pro- 
pose extending this technique to image the relative phase 
profile of an elongated pscudospin- 1 BEC to study the 
spatial dependence of relative phase fluctuations. 
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